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Abstract 

The Low Energy Neutron Detector Array (LENDA) is a neutron time-of-flight 
(TOF) spectrometer developed at the National Superconducting Cyclotron Lab- 
oratory (NSCL) for use in inverse kinematics experiments with rare isotope 
beams. Its design has been motivated by the need to study the spin-isospin 
response of unstable nuclei using (p, n) charge-exchange reactions at intermedi- 
ate energies ( > 100 MeV/u). It can be used, however, for any reaction study 
that involves emission of low energy neutrons (150 keV - 10 MeV). The array 
consists of 24 plastic scintillator bars and is capable of registering the recoiling 
neutron energy and angle with high detection efficiency. The neutron energy is 
determined by the time- of- flight technique, while the position of interaction is 
deduced using the timing and energy information from the two photomultipliers 
of each bar. A simple test setup utilizing radioactive sources has been used to 
characterize the array. Results of test measurements are compared with simu- 
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lations. A neutron energy threshold of < 150 keV, an intrinsic time (position) 
resolution of ~ 400 ps (~ 6 cm) and an efficiency > 20 % for neutrons below 4 
MeV have been obtained. 

Keywords: neutron detector, neutron Time-of-Flight, low energy neutron 
detector, inverse kinematics, (p, n) charge-exchange 



1. Introduction 

A high-efficiency position-sensitive scintillator array for the detection of low 
energy neutrons has been developed at the National Superconducting Cyclotron 
Laboratory (NSCL). The Low Energy Neutron Detector Array (LENDA) was 
built to study charge-exchange (CE) (p,n) reactions. These reactions have been 
successfully used in the past to study the spin-isospin response of stable nuclei 
[TJ [2] . LENDA makes possible the extension of such studies to unstable nuclei 
using the (p,n) CE reaction in inverse kinematics [3 . Moreover, LENDA can 
be utilized in any nuclear reaction study that involves production of low-energy 
neutrons, for example, for proton-transfer ((i,n) reactions in inverse kinematics 
and beta-delayed neutron emission experiments. In such investigations, the 
detection of the slow neutrons with good efficiency as well as energy and angle 
resolution is essential. LENDA addresses this need by providing fast timing 
(~400 ps timing resolution), reasonable position resolution (~ 6 cm ) and a low 
neutron detection threshold of 150 keV. 

For the characterization of the array, a simple test setup utilizing neutrons 
emitted by a ^^^Cf radioactive source and photons from a ^^Na source has been 
used. In this work, the results of the characterization are presented. The setup 
is described in section [2] The electronics and data acquisition system (DAQ) 
is presented in section [3) a brief description of the data reduction process in 
section |4j the efficiency, timing and position resolution of LENDA in sections [5j 
|6] and [7| respectively. The conclusion follows in section |8] 

2. Description of LENDA 

LENDA is an array of 24 neutron detectors. Each detector module is a type 
BC-408 [4 plastic scintillator bar with dimensions of 300 x 45 x 25 mm. A Hama- 
matsu II6410 photomultiplier (PMT) assembly with a photo-electron gain of the 
order of 10^ is used at each end of the bar to detect the scintillation light. In 
each module, to ensure optimum transmission of light to the photocathode, the 
PMTs are directly coupled to the scintillator using optical epoxy. The scintilla- 
tor bar is wrapped with one layer of white nitrocellulose membrane filter paper 
from Advantec 0, a layer of aluminum foil and finally black insulating tape to 
ensure proper light propagation through the bar as well as light-tightness. The 
Advantec film has comparable performance to other high refiectivity scintillator 
wrapping materials available commercially (Gore DRP [6 , 3M Vikuity [7 ) and 
is easier to manipulate while wrapping the LENDA scintillators. 
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Figure 1: One half of the LENDA array (12 detector modules) mounted on its 
stand. In this configuration, optimized for charge-exchange (p,n) experiments 
in inverse kinematics, the bars are arranged at a radius of 1 m from the target 
position. One stand is placed on each side of the beam line and the array covers 
roughly 45 degrees of scattering angle in the laboratory frame. 

Figure [l] shows a picture of one half of the array. Each bar is mounted so that 
the position-sensitive direction of the bar is vertical, and the shortest side of the 
bar is parallel to the neutron fiight path. At a distance of 1 m from a neutron 
source, one array covers a solid angle of 0.16 sr in total for a scattering-angle 
coverage of about 45 degrees. 

By combining the time and pulse-height information from the PMTs, the timing 
of a neutron hit, the corresponding scintillation light output, and the neutron's 
hit position along the longest side of each LENDA bar can be determined. Dur- 
ing standard Time-of-Flight (TOF) operation mode, LENDA has to be used 
with an external detector that provides a time reference signal. In a typical 
experimental setup where LENDA is used in combination with the S800 spec- 
trometer of NSCL [8 , that signal is generated by charged particles hitting a 
diamond detector [9 at the object of the spectrometer. 

3. Data Acquisition system 

A common-stop trigger logic is employed for the time-of-fiight measurement. 
The trigger logic is realized using a Xilinx Virtex-II FPGA chip [10 implemented 
in a VME-based JTEC XLM72V FPGA module [H^. Two 32-channel CAEN 
V792 charge-to-digital converters (QDCs) [12 are used to digitize the charge of 
the PMT signals. The time information for each event is digitized by two 32- 
channel CAEN V775 time-to-digital converters (TDCs). The PMTs are powered 
by ISEG EHS F030n and 8030n high voltage power supplies. Figure [2] shows a 
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schematic diagram of the LEND A DAQ. To achieve the highest possible signal- 
to-noise ratio from the PMTs and thus maintain a low threshold of neutron 
detection, the PMT voltage is set to be as close as possible to the maximum 
value of -2700 Volts permitted by the manufacturer. Some voltage variation (50 
-100 Volts in most cases) between individual PMTs is employed to obtain a 
rough gain-matching of the signals. The anode signal of each PMT is attenuated 
by a factor of 0.3 using a custom-built 16-channel, two-way splitter-attenuator 
module, before it is used as an input to the DAQ electronics. This is done 
to match the signal's pulse-height to the dynamic range of the digitizers and 
discriminators for the neutron energies of interest. For each PMT, one of the 
split anode signals, after being cable-delayed by ^ 100 ns, is sent to the QDC 
to digitize its charge. The other split signal is sent to a Phillips 7106 leading 
edge discriminator to generate a logic signal that indicates the PMT has fired 
(PMT-event signal). The PMT-event signal is sent to the TDC and serves as a 
start-signal for the TOF digitization of that event. It is also used as one of the 
input signals for the FPGA trigger-logic. One more input signal for the logic 
is required to identify a "good" event. That signal is provided by an external 
trigger-detector that is experiment-specific. Five logic signals appear at the 
output of the FPGA as a result of the trigger processing: 

• A LENDA-hit signal, which signifies that at least one LEND A bar is hit. 
This is generated by the coincidence of the PMT-event signals from the 
two ends of each LEND A bar. An OR between those coincidences from 
aU LENDA bars gives the LENDA-hit signal. 

• A QDC gate signal provided by the OR of all PMT-event signals. 

• A FAST CLEAR signal generated by the same OR, with a delay of a few 
hundred ns. This clears the digitizers and aborts the processing of the 
signal when the DAQ is not triggered by a good event. 

• A VETO signal that eliminates the FAST CLEAR signal whenever a co- 
incidence is realized (DAQ trigger) between the external trigger and the 
LENDA-hit signal. This allows "good" data to be digitized successfully. 

• A COMMON STOP signal that causes the TDC to digitize the timing 
information. This is, in terms of timing, the same as the DAQ trigger 
signal. 

The FAST CLEAR scheme that ahows events to be selected after the QDC 
integration has an advantage. The starting time of the QDC integration can be 
set independently of the timing of the trigger signal. This reduces the amount 
of delay required to place the detector signal inside the QDC gate. 

4. Data Reduction 

Neutron TOF, hit position along the longest side of a scintillator bar, and 
light output in the scintillator are obtained by combining the timing of signals 
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Figure 2: A schematic diagram of the DAQ system for LEND A. 



from the PMTs at the top and bottom ends {tr and ts) oi the scintihator and 
the signal amphtudes {qr and qs) recorded in the TDC and QDC modules, 
respectively. Herein, we briefly describe the definition and derivation of these 
parameters. This analysis is similar as described by other authors [13| JA . 
Owing to quenching effects and the propagation time of scintillation light from 
the hit position in the scintillator to the PMTs, the signal timing and amplitude 
are both asymmetric between the top and bottom PMTs. This asymmetry is 
correlated with the hit position. Two position-dependent parameters xt and 
xq are derived from the difference between It and t^, and between qT and qs 
using the relations, 

XT CX Tdiff =tT -tB, (1) 

qr - qs 

XQ (X cogQ = . (2) 

Qt^Qb 

The proportionality between xq and cogQ in equation |2] is a valid approximation 
in the case of LENDA: The length of the scintillator bars (30 cm) is much smaller 
than the light attenuation length of 210 cm for the BC-408 material |4j; then, 
the exponential dependence of qr and qs on hit position can be approximated 
by its Taylor expansion leading to a linear relationship. 

The neutron TOF {t) and light output (/), are derived from an average of It 



5 



o 

^ 1000 



800 



600 



400 



200 




500 1000 1500 2000 2500 3000 3500 4000 
Pulse height (QDC Channels) 



Figure 3: QDC spectrum of neutrons emitted by a ^^^Cf source at a distance of 
Im from the detector. 



and t^, and from qt and qb^ respectively, as 



t = 



I 



(3) 
(4) 



where tTR is the time reference (TR) for TOF provided by the external detector. 
The qt and qb values are converted from QDC channels to light output units 
of electron equivalent energy (MeVge), by assuming a linear relation between 
the QDC values and the light output. This approach has been shown to re- 
produce quite accurately the experimentally measured light output for organic 
scintillators [15 . A typical QDC spectrum for neutrons emitted from a ^^^Cf 
source is shown in figure |3] The relation connecting QDC values to light output 
is calibrated using the 59 keV 7-ray of ^^^Am and the Compton edge of the 
662-keV 7-ray from a ^^^Cs source. 



5. Neutron spectrum and efficiency with Cf 

Time-of-fiight spectra of neutrons from a ^^^Cf fission source were measured 
by LENDA to determine its neutron detection efficiency. The kinetic energy 
of fission neutrons from ^^^Cf ranges approximately from 100 keV to 10 MeV 
according to the IAEA evaluated neutron energy distribution [161 113 • This dis- 
tribution was used to reliably determine the energy dependence of the efficiency 
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Figure 4: A schematic view of the experimental setup used in the measurement 
of a ^^^Cf neutron spectrum. The ^^^Cf fission source was placed at a distance of 
1 m from the LEND A bar. Prompt 7-rays from fission were detected in a liquid 
scintillator. Neutron-7 ray pulse shape discrimination was used to separate them 
from prompt neutrons. The 7-rays identified in the liquid scintillator served as 
the time reference for the neutron TOF. The shadow bar was used to determine 
the contribution from background events due to wall-scattered neutrons. See 
the text for details. 



curve. 

In the experimental setup (shown schematically in figure [4|, the ^^^Cf source 
was placed at a distance of 1 m from the LENDA bars. The TR signal for 
the neutron TOF was established by detecting prompt 7-rays from fission in 
a 2" X 2" cylindrical NE-213 type liquid scintillator manufactured by EL JEN 
(EJ-301) [18 . This detector was placed at a distance of 6 cm from the source. 
Pulse shape neutron-7 discrimination was employed in the TR scintillator for 
selecting only 7-ray events. In all TOF spectra a correction for the "walk" of 
the leading edge discriminator was applied in the off-line analysis. Figure [sja) 
shows a 2-dimensional spectrum of light output versus TOF for 7 rays and neu- 
trons detected in a LENDA bar. The continuous curve in figure (sjc) shows 
the projection of the 2-dimensional spectrum on its horizontal (TOF) axis. In 
figure [sj a), the events included in the sharp peak at t ~ 3 ns are due to 7 — 7 
coincidences between the LENDA and TR detectors. The events that lie on the 
right side of the spectrum at t > 20 ns correspond to 7 — n coincidence events 
(neutrons detected in LENDA). Events appearing between those two regions 
are due to 7 rays from the source, reaching LENDA bars after being scattered 
from the surrounding walls. Any events in the region of t < 3 ns on the left of 
the 7 ray peak come from neutrons leaking through the n-7 discrimination gate 
that were being detected in the liquid scintillator in coincidence with 7 rays in 
LENDA. 
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Figure 5: Response of LEND A to prompt neutrons and 7-rays from ^^^Cf fission 
events shown as a function of the TOF and light output in the scintillator bar 
with (a) and without (b) a shadow bar. The projection onto the TOF axis of 
(a) and (b) is shown in (c). 

Since increasing TOF corresponds to lower neutron kinetic energies, the max- 
imum detected light output in a plastic scintillator like LENDA is varying as 
a function of TOF. The characteristic light output curve corresponding to full 
energy deposition by neutron events, becomes practically zero for TOF values 
of t ^ 150 ns. The TOF for this minimum amount of light output corresponds 
to the smallest kinetic energy that can be detected by LENDA, about 150 keV. 
It is defined by the lowest threshold setting of the DAQ discriminator module 
and the maximum gain of the PMTs. 

The events lying above the light output curve are background events, which are 
mostly due to neutrons and gammas indirectly reaching LENDA after scattering 
by the surrounding walls or objects. For these events, TOF is not correlated 
with the deposited energy. In the offline analysis, the contribution from this 
background can be reduced by eliminating the events in the region above the 
light output curve (light-output cut). 

The remaining background contribution was directly measured by inserting a 
copper shadow bar to block neutrons directly coming from the source to LENDA 
as shown in figure [4j The dimensions of the shadow bar were selected so that 
they match the size of the scintillator bar. It was thick enough (30 cm) to at- 
tenuate neutrons with kinetic energy up to 10 MeV, by a factor larger than 600. 
The effect of the shadow bar blocking part of the background neutron ffux was 
calculated assuming an isotropic background neutron distribution. The shadow 
bar concealed only 5% of the total solid angle subtended by LENDA and the 
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systematic error induced in the efficiency calculation was estimated to be of 
the order of 1%. The live time of the measurement was used to normalize the 
background runs to those without the neutron blocker. The spectrum measured 
with the shadow bar is shown in figure jsjb) , where most of the neutron events 
are due to scattering from walls. The dashed curve in figure [sjc) corresponds 
to the projection of this 2-dimensional spectrum on the TOF axis. 
To determine the absolute value of the efficiency, the ^^^Cf neutron yield distri- 
bution of Refs. [I6j [17] was normalized using a yield reference (YR) detector of 
the same type and size as the TR one. The YR detector was placed at a distance 
of 40 cm from the source. Its efficiency was simulated using the computer code 
NEFF7 [19] and found to vary between 30 and 40% for neutrons with energies 
in the region of 1 to 3 MeV. According to previous studies using fission fragment 
detectors [20], the results of efficiency simulations using the NEFF7 code are 
known to be reliable within an error of 5% for this type of detector. 
After subtracting the contribution from the background, the number of detected 
neutrons, N{En)^ as a function of neutron energy, En^ are compared to the neu- 
tron yield from the ^^^Cf source, Y{En) using the relation, 

N{En) = Y{En) X e{En) X Al], (5) 

where A^l is the solid angle of the LENDA bar precisely calculated from the 
alignment information of the detector and source, Y{En) is the standard shape 
of neutron energy distribution taken from Refs. j6l[T7] and normalized using the 
YR detector, and the e{En) is the neutron-detection efficiency of the LENDA 
bar as a function of E^- The obtained efficiency curves are shown with dots 
in Fig. [6) where three threshold levels of 30, 50, and 60 keVge on the light 
output were applied in the offiine analysis. The results of the simulation of 
the setup using Monte Carlo N-Particle Transport Code (MCNPX) [21 are 
shown for comparison. The thresholds for the recoiling proton energy in the 
simulations were set to be 200, 300, and 400 keV, respectively, to match the 
light-output threshold in the experimental data. The simulated distributions 
reproduce fairly well the trend and magnitude of the measured efficiency for 
the three values of the threshold. Any small deviations are probably due to 
ambiguities in background subtraction. 

6. Timing resolution 

The timing resolution of LENDA was determined using the two correlated 
511 keV photons from the positron annihilation of ^^Na. Due to compton- 
scattering (the dominant mechanism of photon detection in BC-408 at these 
energies) photons of 511 keV produce a continuum of low-energy events which 
are equivalent (in terms of pulse height) to the proton recoils produced in the 
scintillator by 2.3 MeV neutrons. The timing response of the BC-408 plastic 
scintillator to gammas and neutrons is therefore the same (same pulse shape) 
with the pulse- height of the compton edge of 511 kev photons corresponding to 
a neutron energy of 2.3 MeV. 
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Figure 6: Experimental neutron detection efficiencies of a LEND A bar (dots) as 
a function of neutron energy in comparison with the curve simulated by MCNPX 
(curves). The error bars (shown where they are larger than the points) indicate 
statistical uncertainties. See the text for details. 




Figure 7: Schematic depiction of the setup used to determine the timing reso- 
lution of LENDA bars. Photons from a ^^Na positron source were used. In this 
setup, one of the LENDA bars was functioning as the trigger detector providing 
the time reference signal for the DAQ. 



For the timing resolution measurements, two LENDA bars were placed oppo- 
site to each other so that the source would be at a distance of 1 m from one of 
the bars, as shown in figure (tI One of the LENDA bars (closest to the source) 
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served as a trigger detector while the other registered the photon TOF. In the 
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Figure 8: Spectrum representing the timing resolution obtained with LEND A 
using a ^^Na positron emitting source in the arrangement shown in figure [?[ The 
peak corresponds to 7-7 photon coincidence events from positron annihilation 
that are detected by the two LENDA bars. The width of the peak corresponds 
to a timing resolution of ^420 ps. 

time spectrum of figure [8] a prominent peak corresponding to 7-7 coincidences 
is observed. The width of the photon peak in figure [S] is roughly 600 ps FWHM 
which corresponds to an intrinsic resolution of 600 ps/a/2 ^ 420 ps for each 
bar. This timing resolution is nearly constant for light outputs above 200 keVge 
and degrades fast for smaller light outputs. As shown in figure |9) this timing 
resolution corresponds to an almost constant 5-6% /S.E/E resolution for neu- 
tron energies below 4 MeV when detected by LENDA at a distance of 1 m. 
The degradation of timing resolution for light outputs below 200 keVge causes 
to the energy resolution only a small deviation from linearity for low energy 
neutrons (below roughly 1.3 MeV). In an case this degradation of the timing 
resolution does not significantly affect the energy resolution of the detector at 
the corresponding energies. 

7. Position Resolution 

Using the TDC and QDC information from each PMT as described in sec- 
tion |4j the position of an event registered by LENDA can be determined. The 
PMT closer to the event registers the scintillation light sooner and produces 
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Figure 9: Plot of the energy resolution achievable by LEND A for neutron ener- 
gies in the range of - 4 MeV, based on measured timing resolution for photons 
of 511 keV. 



a larger signal than the PMT further away. The position resolution therefore 
depends on the distance traveled by the scintillation light inside the detector 
as well as the attenuation caused to the scintillation light while traversing the 
scintillator material and interacting with the reflective coating of the scintilla- 
tor. The amount of light produced by an event is also a significant factor and 
the position resolution is expected to be better for higher deposited energies in 
the scintillator. 

The position resolution of LENDA bars was determined experimentally using 
radioactive sources. In earlier investigations using a prototype detector and a 
^^^Cf source, a clear correlation was observed between the quantities x^ and 
xq described in equations [l] and [5] [3 . In figure 10, the 2-dimensional spectra 
of those two parameters are presented for a source placed at different locations 
along a detector. The clear shift of the locus of events is directly correlated to 
the change in location of the source. 

For the position resolution measurements, a bore hole in a large lead block was 
used to collimate 7 rays from the ^^Na source. A schematic of the setup is shown 
in figure 11 The thickness of the lead (~10 cm) was enough to block the 511 keV 
7 rays emitted from ^^Na. The diameter of the collimator hole was 1.6 cm. The 
collimator was placed at the minimum distance possible from the detector - 5 
cm - to minimize the diameter of the irradiation area. The resulting irradiation 
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Figure 10: Events from a source attached on the long side of a LEND A bar 
generate a locus in the 2-D spectrum of time difference versus charge asymmetry. 
In this figure the two axes are labeled Xt and Xq respectively and calibrated 
in units of distance (m) . The locus of events changes with source location along 
the scintillator bar, revealing the position sensitivity of the detector. In each 
figure the distance of the source from the middle point (center) of the bar is 
displayed as well. 



spot of 2.4 cm total diameter, was small compared to the position resolution 
expected by such detectors (several cm FWHM). In this setup, LENDA was 
moved vertically in front of the collimated source for the measurements. The 
position of the source center along the bar was read from a scale fixed on the 
side of the detector. 
For five different locations of the collimated source, Tdiff and cogQ (defined 
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Figure 11: Picture of the setup used for position calibration 
above in equations [l] and [2] respectively) between the 2 PMTs of each bar were 



recorded. Typical spectra of these quantities are presented in figure 12 The 



deviation from a Gaussian shape for the distributions in this figure comes from 
photons scattering on the support structure. In the analysis, these scattering 
tails where not considered as valid events and a Gaussian fit was used to deter- 
mine the centroid values of Xq and X^ distributions. The results correlating 
the position of the source with respect to each one of the two parameters are 
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Figure 12: Plot of the T^iff (top figure) and cogQ (bottom figure) distributions 
that correspond to one location of the collimated source along a LENDA bar. 



presented in figures 13(a) and |13(b)] respectively. No constraint was set on the 
deposited energy of the events and therefore these results correspond to an av- 
erage position resolution. A linear trend is prominent in both graphs. The 
small but negligible deviation from linearity observed in the plot of figure |13(b)| 
is attributed to the presence of edge effects that affect the pulse-height of the 
signal for events very close to the two extremes of the scintillator. To better 
understand if the transport of light in the scintillator bars can account for this 
effect, simulations were performed using the code GUIDE? [22]. The code was 
modified to correctly reproduce the transport of light in the area between the 
scintillator material and its wrapping and to include as a parameter the reflec- 
tivity of the wrapping material [23 . In figure 14, the comparison between the 
simulation results and the experimental data of figure p!3(b)| are presented. The 
results of the light transport simulations reproduced fairly well the position de- 
pendence observed experimentally. An average position resolution of 6.3±0.2 
cm was achieved using the timing information (equation 3) while a slightly worse 
resolution of 7.2±0.3 cm was achieved using the QDC charge signal. 



8. Conclusion 

LENDA, an array of plastic scintillators capable of detecting low energy 
neutrons with high efliciency and angular resolution has been built and char- 
acterized at NSCL. In this work the results of the characterization of the array 
using radioactive sources, are presented. The efliciency as well as the timing 
and position resolutions of LENDA have been determined experimentally. The 
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Figure 13: Position dependence of (a) the time difference and (b) of the charge 
center of gravity between the two PMTs of a LEND A bar. The straight hnes 
in the graphs correspond to a hnear fit of the data. The error in the horizontal 
axis corresponds to the diameter of the aperture used to cohimate the source. 

measurements have been supported by Monte Carlo simulations and general 
agreement was observed between data and calculations. 

A low neutron energy threshold of ~ 150 keV has been demonstrated for LEND A. 
A timing resolution of ~400 ps and a position resolution of ~6 cm have been 
obtained. An efficiency greater than 20% for neutrons below 4 MeV has been 
measured in the test setup. The results of this work confirm that LENDA is an 
adequate tool for the study of nuclear reactions in inverse kinematics, especially 
when low-energy neutron detection is required. Recently, the LENDA array 
has been used at NSCL to study for the first time the spin-isospin response 
of the radioactive nucleus ^^Ni [24J. This was the first experiment with a new 
measurement technique of (p,n) charge-exchange reactions in inverse kinematics 
at intermediate energies on unstable isotopes. The results of this investigation 
proved to be important for improving the description of electron-capture rates 
on nuclei in the iron region, and for modeling the late evolution of core-collapse 
and thermonuclear supernovae. 
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Figure 14: Comparison between simulations of the light transport in LENDA 



with the code GUIDE? (red squares) and the experimental data of figure 13b 



(black dots). The error in the experimental data due to the aperture size is 
not included explicitly in this graph but the increased size of the experimental 
points is used to point out the existence of experimental uncertainty. The dashed 
straight line is a fit excluding the 2 outermost experimental points, used to guide 
the eye to the deviation from linearity at the edges of the detector. 
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